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o. 

, Star formation rates in the centers of disk galaxies often vastly exceed those 

at larger radii, whether measured by the surface density of star formation Eg F R, 
Ck by the star formation rate per unit gas mass, Esfr/E, or even by total out- 

put. In this paper, we investigate the idea that central starbursts are self- 
regulated systems, in which the momentum flux injected to the interstellar 
medium (ISM) by star formation balances the gravitational force confining the 
ISM gas in the disk. For most starbursts, supernovae are the largest contrib- 
utor to the momentum flux, and turbulence provides the main pressure sup- 
port for the predominantly-molecular ISM. If the momentum feedback per stel- 
lar mass formed is p*/m* ~ 3000 km s -1 , the predicted star formation rate is 
E SF r ~ 27rGE 2 m*/p* ~ 0.1 M kpc -2 yr -1 (E/100 M pc~ 2 ) 2 in regions where 
gas dominates the vertical gravity. We compare this prediction with numer- 
ical simulations of vertically-resolved disks that model star formation includ- 
ing feedback, finding good agreement for gas surface densities in the range 
E ~ 10 2 — 10 3 M pc -2 . We also compare to a compilation of star forma- 
tion rates and gas contents from local and high-redshift galaxies (both merg- 
ers and normal galaxies), finding good agreement provided that the conver- 
sion factor Xco from integrated CO emission to H 2 surface density decreases 
weakly as E and Esfr increase. Star formation rates in dense, turbulent gas are 
also expected to depend on the gravitational free-fall time at the correspond- 
ing mean ISM density po; if the star formation efficiency per free-fall time is 
£ff(po) ~ 0.01, the turbulent velocity dispersion driven by feedback is expected 
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to be v z = 0.4£ff(po)p*/^* ~ lOkms -1 , relatively independent of E or Esfr- 
Turbulence-regulated starbursts (controlled by kinetic momentum feedback) are 
part of the larger scheme of self-regulation; primarily-atomic low-E outer disks 
may have star formation regulated by UV heating feedback, whereas regions 
at extremely high E may be regulated by feedback of stellar radiation that is 
reprocessed into trapped IR. 

Subject headings: galaxies: starburst - galaxies: ISM - ISM: kinematics and 
dynamics - ISM: star formation - turbulence 



1. Introduction 

Averaged over large scales in disk galaxies, the relationship between the mean surface 
density of star formation, E SFR , and the mean gaseous surface density, E, is observed to 
be superlinear, bot h in the local Universe and at higher redshift (e.g. iKennicuttl Il998bl ; 
Genzel et al.ll2010al ). This behavior reflects the increased efficiency (or shorter timescale) of 
star formation under conditions of higher mean gas density, which is correlated with higher 
E. Starburst regions within the central ~ kpc of galaxies, commonly ob served as dust- 
enshrouded LIRGs and ULIRGs in which the ISM is predominantly molecu lar (ISanders & Mirabel 
19961 : iGenzel et al.lll998t lKennicuttlll998al : ISolomon fc Vanden Boutll2005l ). represent the ex- 



treme of this behavior, with E ~ 10 2 — 10 4 M Q pc -2 and the star formation (or gas conver- 
sion) timescale tsF, g as = E/Esfr a factor of 10 or more below tsF, g as in lower- E regions of 
galaxies. Although earlier observations based on global averages suggested a simple power- 
law relationship between E and Eg F R, these global relations have considerable scatter about 
the mean. More recently, evidence has emerged that the star formation regime within galac- 
tic centers at very high E ^ 100 M pc~ 2 is distinct from the star formation regime that 
prevails in the main disks of galaxies at lower E £ 100 M pc~ 2 . Within the main disk, an- 
other change in the star formation regime appears to occur from "mid-disk" to "outer-disk." 
Typically, both galactic center and mid-disk regions are primarily molecular, whereas outer 
disks are primarily atomic. Figure ([T]) provides a schematic, dividing the disk into different 
star-forming regimes. 

In main-disk regions where E 100 M pc~ 2 , high- resolution observations have found 
that EgpR is proportional to S mn i , the mean surface density of mol ecular gas, with a timescale 
^SF,moi = E mol /E SFR « 2 Gyr faigiel et al.ll2008l : felanc et al.l2009| ). The value E ~ 100 M pc~ 2 



that defines the upper limit for "main disks" is characteristic of the surface density Eqmc of 



(Blitz et al. 


2007; 


Sheth et al. 


2008; 


Bolatto et al. 


2008) 



The relation E SFR oc E mol for 
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Fig. 1. — Schematic indicating the geography of galactic disk star formation; ranges for 
gaseous surface density £ (in M Q pc -2 ) in the various regions are approximate. Regula tion 
of star formation in the main disk has recently been discussed by lOstriker et al.l (120101 ) . In 
this paper, we analyze self-regulation of star formation in the turbulence-dominated region 
within the galactic center, for starburst disks. 



£ ^, S GMC is consistent with the idea that spatially-isolated GMCs have an approximately 
constant specific star formation rate; this rate, and the value of Somc , may be controlled by 



internal feedback processes within GMCs (e.g. iKrumholz et al.ll2006l ). al though this is no t 
yet completely understood. Over a typical GMC lifetime of ~ 20 Myr ( iBlitz et al.l 120071 ). 
~ 1% of the gas must be converted to stars in order to yield tsF,moi ~ 2 Gyr within main-disk 
regions. Although there are uncertainties arising from limited resolution and calibration of 
the CO-to-H 2 conversion (see Section 4), the relationship between S mo i and Esfr in galactic 
center regions with X ^ 100 M pc -2 appears considerably steeper than linear. Presum- 
ably, this is because the molecular gas in galactic center regions with £ > 100 M pc~ 2 is 
at higher mean (volume) density than the molecular gas in main-disk GMCs having "stan- 
dard" surface density Egmc ~ 100 M pc -2 . 

Molecular gas dominates mid-disk regions, but is a small fraction of the total in outer 



2004. 


2006: 


Lerov et al. 


2008). 



explain the star-forming behavior and balance of phases in the main regions of disks at 
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E 100 M Q pc -2 , based on the requirement that the diffuse H I gas must be in both ther- 
mal and dynamical equilibrium. In this theory, star-forming gravitationally-bound clouds are 
assumed to have internal pressures much larger than their surroundings, consistent with local 
GMC observations. Thus, they are effectively isolated, and over their lifetimes are assumed 
to maintain a mean internal specific star formation rate consistent with the empirical main- 
disk value tsF.moi ~ 2 Gyr. In equilibrium, individual gravitationally-bound clouds form and 
disperse at equal rates, with neutral ISM gas cycling between the diffuse and gravitationally- 
bound components. In outer-disk regions (typically E £ 10 M pc -2 , although this varies 
with p*), the equilibrium fraction of gas in star- forming bound clouds is just enough that 
the heating provided by the resulting stellar UV balances the cooling of the surrounding 
diffuse H I; here, Esfr oc EyZ/v In mid-disk regions (typically 10 £ E £ 100 M pc~ 2 ), 
the interstellar medium becomes predominantly molecular (by mass) because the heating 
provided by star formation would be insufficient to balance cooling in high-density warm 
atomic gas ( which is compr e ssed by the vertical disk gravity); here Eg FR oc E mo i oc E. 
The model of lOstriker et al.l (120101 ) assumes that most of the volume is filled by non-star- 
forming atomic gas, and thus is inapplicable to the starburst regime in galactic centers where 
E 100 M Q pc -2 and gas is mostly molecular. 

Where E is very high, the conditions throughout the ISM resemble those in the interior 
of locally-observed GMCs, but are even more extreme in terms of the mean densities, the 
abundance of very high de nsity gas, and the ratio of the turbulent velocity dispersion Vt m h to 



the thermal speed y^u (e.g. IScoville et al.lll99ll ; lDownes et al.lll993l ; lDownes fc Solomonlll998 



Gao fc Solomon!l2004l ). Supersonic turb ulence is known to be a dominant process controllin g 
star formation regulation within GMCs (IMac Low fc Klessenll2004j ; iMcKee fc Ostrikerll20071 ). 
so it is expected to be important in starburst regions as well. Main-disk GMCs are relatively 
isolated entities, in the sense that their internal dynamical timescales are short compared to 
the orbital, epicyclic and vertical oscillation times in the galactic potential (fi -1 , k" 1 , and 



respectively). In galactic center regions, however, the gravitational potential is deeper, 
such that Q, k, and v are all substantially higher, with timescales ^ 10 Myr. Thus, the 
galactic environment may also significantly impact the development of star formation. 

Turbulence in the ISM both inhibits and encourages star formation. For regions large 
compared to the energy (and momentum) injection scale, turbulence provides an effective 
pressure that opposes gravity. In the ISM, this helps to regulate the mean density av- 
eraged over the disk thickness p, and therefore the large-scale self-gravitational timescale 
(tgrav) oc p -1 / 2 oc (/f/E) 1 / 2 . Here, the gas disk semi-thickness H varies with turbulent 
velocity dispersion as H oc a z or cr 2 depending on whether stellar or gas gravity domi- 
nates the vertical potential gradients. In either case, larger turbulent velocities raise H 
and (tgrav) for the ISM disk, while larger stellar and gas densities lower H and (t gTav ). On 
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scales smaller than the turbulent injection, supersonic compression creates local overdense 
regions pi oca i ^ P that potentially can collapse more rapidly than the larger-scale regions 
containing them since £ graV) i oca i/(t grav ) oc (pi oca i/ 'p) 1 ^ 2 ^ 1- However, turbulent rarefactions 
and shear also destroy overdensities. Because t graVi i cai depends on density, turbulence will 
disperse the moderate-density structures before self-gravity can concentrate them further, 
while allowing the highest- density regions (amounting to a small fraction of the mass) to col- 
lapse and form stars. Theoretical arguments suggest that the fraction of mass in turbulent 
systems that collapses to f orm stars will be relatively insensitive to the velocity dispersion 



( IKrumholz &: McKed l2005t see also Section 5), so that the overall star formation rate is 
primarily governed by the gravitational time on large scales, (t grav ). Since massive stars 
energize the ISM, raising the velocity dispersion and (t gTav ), the associated feedback loop 
may allow star formation rates to be self-regulated. 

In this paper, we investigate the process of self-regulation via turbulent driving, for 
application to understanding what controls star formation rates in starburst disks. For 
molecule-dominated regions, cooling times are short, so it is the injection of momentum, 
rather than energy, that is essential in defining the dynamical state of the disk. Thus, we 
begin in Section 2 by considering the implications of maintaining force balance in the ver- 
tical direction, which imposes a requirement on the input momentum flux associated with 
star formation feedback. We show that for disks where supernovae dominate the feedback, 
the star formation rate is expected to vary as Esfr oc £ 2 . At very high surface density in 
optically-thick disks, radiation pressure can exceed the turbulent pressure driven by super- 
novae, which leads to S SFR oc S. In Section 3, we compare our prediction for S SFR vs. S 
in turbulence-dominated disks to the results of numerical simulations, for which feedback is 
explicitly implemented in local vertically-resolved models of self- gravitating, rotating disks. 
Section 4 compares the prediction for Esfr vs. X to observations of galaxies in the regime 
where the turbulence-dominated model is expected to apply. Both the numerical simula- 
tions and the observations are consistent with the simple analytic model. In Section 5, we 
connect to models relating star formation rates to (t grav ), and discuss how the mean inter- 
nal properties (including velocity dispersion and disk thickness) of turbulence-dominated, 
feedback-regulated disks are expected to depend on E and the parameters that characterize 
the star formation feedback. We conclude in Section 6 with a brief summary and discussion. 
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2. A simple model for momentum-controlled, self-regulated disk star 

formation 

Consider a region of a galactic disk with gas surface density E, residing within a bulge of 
uniform density p&, such that the gas-disk gravitational potential is $ ff and the bulge potential 
is We assume that the gas disk's internal structure is similar to that of a GMC, in that it is 
highly inhomogeneous (due to supersonic turbulence) but most of the gas concentrations are 
transient. The disk is thus taken to consist primarily of "diffuse" gas, in the sense that only 
a small fraction of the mass at any time lies in bound clumps having gravitational potential 
large compared to the mean midplane value in the disk. With an effectively plane-parallel 
distribution of mean gas density, the mean gravitational potential depends only on the 
distance from t he midplane z. The time-averaged, horizontally-averaged vertical momentum 
equation (e.g. iBoulares &: Cox 1990 ; Piontek Sz Ostriker 2007 ; Koyama &; Ostriker 2009b ; 
Ostriker et all bold ) may then be vertically integrated, yielding the result that the total 
momentum flux through the disk must be equal to the total vertical weight of the overlying 
gas. 

At the midplane, the weight includes a self-gravitational term, 



dz 



'-dz 



8ttG 



dz 



-dz 



tiGT? 



where we have used p = (AitG) 1 d 2 § g /dz 2 and {d^g/dzloo 
arising from the external bulge potentials, 



27rGS for a slab; and a term 



p— — dz 

dz 



AirGpb 



pzdz = 



3p 



(2) 



Here, p(z) is the mean (area-weighted) density of the gas at height z, and p is the midplane 
value. The value of the numerical coefficient Q is insensitive to the exact vertical distribu- 
tion of the gas, equaling 1/tt for a Gaussian (when ${, dominates) and ln(2)/2 for a sech 2 
distribution (when <3> 9 dominates); Q ~ 0.33 within 5% for these distributions. Note that 
the self-gravitational weight is independent of the midplane gas density, whereas the weight 
in the external bulge potential is proportional to p&/ p . 

The difference in momentum flux between the disk midplane and the height z max at 
which p — > includes thermal and turbulent terms, P t h = Po( v th) an d -Pturb = Po(vl), where 



1 Note that the expre ssion ( [2| for a stellar bulge differs from the corresponding equation for a stellar disk 
lOstriker et aD |2010l ) by a factor 1/3. 
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(v 2 ) is equal to 1/3 of the total turbulent velocity variance v 2 urh if it is isotropic. The 
thermal and turbulent terms can be combined as a single midplane kinetic pressure po^ 2 , 
where a z is the observed velocity dispersion for a face-on disk. Note that P th corresponds 
to a volume-weighted mean pressure, and (v 2 h ) corresponds to a mass-weighted average of 
thermal speeds of different components: if M tota i = J2i Mi and V to tai = Yli V% for mass and 
volume with Pl = Mi/Vi = P^;J<v then P th = £. V t P t h,i/V total = E i M l v 2 hi /V tc 



total 



PoJ2i(Mi/M tota i)vt hi = p (v 2 h ). For niolecule-dominated regions, the RMS thermal speed 
(vth) 1 ^ 2 ~ 1 km s _1 if T £ 100 K, whereas observed velocity dispersions are at least several 
km s -1 so that (v 2 h ) <C {v 2 z ) w <j\. For convenience, we shall use v z to denote the RMS value 
(v 2 z ) 1/2 - 



The momentum flux difference also contains magnetic terms, AP 1 



mag 



A(P 



B 2 )/8tt, where A denotes the difference between midplane values and values at the sur- 
face of the neutral-gas disk. Numerical s i mulations of highly-s upersonic magnetohydrody- 



namic turbulence (e.g. IStone et al.lll998t iMac Low et al.lll998l ) have shown that the mag- 
netic energy is typically smaller than the turbulent energy by a factor ~ 2 — 3, so that 
B midplane / ft 71 ) ~ Po{ v z)') this * s a ^ so consistent with observations ( iFerrierd 120011 ) . Since the 
magnetic scale height L B is typically large compared to the semi-thickness H = E/p of 



the neutral-gas layer in galactic disks (jTullmann et al.ll2000l ; iHeesen et al.ll2009l ). the ratio 
AP mag /P ma g ~ H/Lb <C 1, so that the resulting magnetic contribution to the momen- 
tum flux difference will be small ( £ [H / Lb\po(v 2 )) . Cosmic rays and radiation may also 
contribute to the momentum flux difference, with respective terms AP cr and AP ra( }. In- 
cluding all terms, the total momentum flux corresponds to an effective midplane pressure 

Peff = -Pturb + -Pth + AP mag + AP cr + AP ra( j. 

By equating P e ff with the total weight (the sum of equations [1] and [2]) , we obtain 



P 



eff 



-fturb 4 

ttGE 2 



P th + AP mag + AP cr + AP rad = p a 2 z (l + K) 



;i+x) 



Here, 



n 



AP mag + AP cr + AP, 



rad 



(3) 



(4) 



characterizes the contribution to vertical support from non-kinetic compared to kinetic terms. 
The term \ = 4CdPfe/(3po) represents the contribution to the gas disk's weight from the 
external bulge potential compared to the gas disk's self-gravitational weight. Since the 
midplane gas density is given by p = 4£ d p 6 /(3x), the vertical equilibrium equation (J3J) is a 
quadratic in \ that can be solved to obtain 

2C(1 + K) 



X 



[1 + 4C(1 + K)} 



1/2 : 



(5) 
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where C = 8Qpb(J 2 / \3iiGT, 2 ) depends only on large-scale properties of the system. Note that 
C 1 / 2 is approximately equal to the ratio of scale heights for a disk-only potential compared 
to a bulge-only potential, so that C <C 1 if the gas disk dominates vertical gravity and C ^> 1 
if the bulge dominates vertical gravity. 

Equation ([5]) can also be rearranged to yield C = x(l+x)/(l+7£); when TZ, C -C 1, x ^ 
C. Using p b = 3fi 2 /(47rG), C = 2QVV 2 « 0-7VV 2 in t erms of the quantity W = (<7^)/(ttGE). 



The parameter W is the analog of the iToomrd (119641 ) parameter Q = kurJ (7tGS) for vertical 
rather than horizontal oscillation frequency and velocity dispersion; i.e. W = Qa z j{2a^) for 
epicyclic frequency k = 2Q. Thus, if Q ~ 1 — 2, TZ -C 1, and (j z /gr £ 1, then C £ 0.2 — 
0.7 and x ~ 0.2 — 0.5 so that self-gravity dominates the external (stellar bulge) gravity in 
controlling the vertical pressure, and p £ p b . Alternatively, 



Q 





r 




1/2 






1 + ^ . 





(6) 



so that for TZ 1 and x ~ 1) Q ~ Sx 1 ^ 20 "/?/ "^- Also, note that the bulge potential generally 
dominates orbital motion (even if the gas potential dominates vertical motion), with a factor 
~ (R/H)C for the stellar vs. gas contributions to Q 2 . 

If star formation is self-regulated via momentum inputs, we expect the turbulent portion 
of the vertical momentum flux through the disk, po(v 2 ), to be comparable to the total 
vertical momentum per unit time per unit area that is injected in the gas by feedback from 
star formation. For isotropic momentum injected per massive star near the midplane, 
averaging over spherical shells yields a vertical component injected on each side of the disk 
of p*/4. We shall therefore adopt the relation 



■Pturb = Po(v 2 z ) = f P 



p* 



4m* 



J SFR) 



(7) 



where m* is the total mass in stars formed per massive star. The ratio p^/m^ is simply 
the mean radial momentum injected in the ISM per unit mass of stars formed. We have 
introduced a factor f p to parameterize the dependence on the details of turbulent injection 
and dissipation. The value of f p is expected to vary between 1 (for strong dissipation) and 
2 (for weak dissipation; this is the limit for a stream of particles that is injected vertically 
and conserves energy as the particles fall back to the midplane). 

The kinetic momentum injected in the ISM per massive star can have contributions 
from a number of sources, inc l uding expanding H II regi ons, stellar winds, and supernovae 
(e.g. iNorman &: Ferraralll996l ; iMac Low &: Klessenl 12004 ) . If star clusters are born in suffi- 



ciently optically-thick regions with deep gravitational potential wells, ra diation pressure also 



becomes important in accelerating the residual gas to high velocities (IMurray et al.l 12010 
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Krumholz fc Matznerl 120091 ). iMatznerl (120021 ) estimates the momentum injected by an ex- 
panding H II region for a given source of ionizing luminosity and ambient density within a 
GMC, finding a ratio of momentum-to- mass p^/m * ~ 20 — 300 km s _1 for clouds of mass 
~ 10 5 — 10 7 M . Using the results of IMatznerl (120021 ). it can be shown that p^/m* oc 



^GMC-^GMC e GMC 101 ^gmcj -Rgmc, and £gmc the mass, radius, and integrated star- forming 
efficiency of a GMC, so that GMCs with escape speeds (GM GMC / Rqmc) 1 ^ > 9 km s _1 
or Eqmc > 0.01 would have p*/m* £ 200 km s _1 . Radiation pressure is most important 
when gas is concentrated near a central stellar cluster and the efficiency is high (see Ap- 
pendix A), so that the value of p*/m* will be comparable to the escape spe ed from the (su- 
per) cluster that forms - up to ~ 100 km s" 1 for the most extreme clusters (IMcCrady et al. 
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20031 : iMaraston et all 12004 ; lOverzier et all 120091 ). 



Although radiation and ionized-gas momentum inputs are important for destroying indi- 
vidual gravitationally-bound GMCs (thereby limiting star formation), supernovae a re likely 
the m ost important feedback mechanism for driving turbulence in the ISM as a whole (ISpitzer 
19781 ). For supernova energy of 10 51 erg and ambient density 10 3 — 1 cm" 3 , the momen tum 
injection per event is m 2 — 5 x 10 5 M Q km s -1 (IChevalierl 1 1974 ; ICioffi et al.lll988l ). If 



multiple supernovae combine to drive a single expanding shell, stellar winds contribute as 
well, but this increases the total energy (and moment um, which is approxi mately l inear in 
the in put energy) per unit mass by less than ~ 20% ( Leitherer et al.lll999h . For a iKroupa 
(120011 ) IMF, the total mass in stars per high mass star (M > 8 M©) is m* = 100 M , so we 
shall adopt p^/m* = 3000 km s _1 from supernovae as our fiducial numerical value. 

Under the assumption that cosmic rays are accelerated by supernova blast waves and 
diffuse vertically out of the disk, the momentum flux contribution from cosmic rays will be 
proportional to that from supernovae, with 



cr , 

L cr 2v A m* ' 



J SFR- 



Here, ( cr is the efficiency of cosmic ray production per supernova of energy -E SN , and we have 
assumed that the cosmic ray fluid streams at speed va - The mean free path of cosmic rays, 
L cr , is likely compa rable the magnetic scale height Lb (jYan fc Lazarianl 120081 ). Taking ( cr ~ 
0.1 jRevnoldsll20"08l ). v A ~ 10 km s" 1 , and £ SN = 10 51 erg, {E SN ( cr )/(2v A m*) = 2500 km s _1 ; 
this is comparable to the kinetic input p*/m* from supernovae. Since the neutral-disk thick- 
ness is small compared to the magnetic scale height, however, H/L CT ~ H/Lb 1 and 
the vertical support of the neutral disk from cosmic rays is small compared to that from 
turbulence given in equation (j7j), AP cr /(p (t>^)) ~ H/L cr 1, similar to the situation for 
magnetic support AP mag /(p (^)) ~ H/L B <^ 1. 



Radiation pressure is important to vertical support of the disk if it is optically thick 



- 10 - 



to reprocessed infrared radiation ({Thompson et al.ll2005l ). Assuming uniformly-distributed 
sources at the disk midplane, the radiation pressure term is 

AP rad = = ^f^S SFR . (9) 

Here, k ir is the mean opacity, F rad is the total luminosity per unit area produced by the disk 
(half is radiated from each si de), c is the speed of light, and based on a standard Starburst99 
model (ILeitherer et alJll99flh . e* = Lj(c 2 M * \ ta 6 2 x 10" 4 is the mass-to-radiation energy 
conversion efficiency by stars for a iKroupal (120011 ) IMF in steady state (applicable if the 
starburst duration is ^ 10 7 yr). The radiation pressure term becomes comparable to the 
kinetic term driven by supernovae (eq. [7]) if e*CKiRX! approaches p*/m*. 

Retaining just the turbulence and radiation pressure terms in equation (jHJ), and using 
equations ((Tj) and ©, the surface density of star formation is given by 

2vr (1 + x ) m*G£ 2 



Here, 



fpP* 



J SFR 



16/ P 



fp 1 + r/r, 



V* 



3000 km s- 



6.2 x 10" 4 



and the optical depth through the disk is 



r = k ir £ = 0.21 



10 cm 2 g" 



100 M pc 



-2 



(10) 



(11) 



12) 



Trapped radiation begins to affect the star formation rate when t/t* £ 1, for very high 
gaseous surface densities^] Note that in equation (jJJ), 1Z — )■ t/t* if radiation dominates over 
the magnetic and cosmic-ray terms and turbulence dominates the thermal term. 

In the turbulence-dominated regime (t/t* 1) for self-regulated disks, the surface 
density of star formation is 



^SFR,turb 



2tt (1 + x) m*GY? 



u 



p* 



0.092 M Q kpc" 2 yr" 



1 + x) ( P*/m* 



3000 km s" 



100 M pc 



-2 



• (13) 



2 Streaming radiation, as well as diffusing radiation, injects momentum to the disk when it is absorbed. 
The maximum momentum flux that can be injected to each side of the disk by streaming stellar radia- 
tion (including ionizing and far-UV radiation before it is reprocessed to IR) is the input value, _F' ra( j/(2c). 
The maximum effect of streaming radiation on self-regulating disk star formation is obtained by replacing 
kir£/2 = t/2 by 1 in equation ©, so that t/t* — > 2/r* = 2e*m :t c/ (f p p r ) in equation ([TU]) . Using the 
fiducial value of t* given in equation (jlip , the maximum effect of momentum input from streaming stellar 
radiation is ~ 10% compared to that from supernovae. 
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As noted above, disks that are marginally gravitationally unstable (Q £ 2) have \ £j 0.5. 
Thus, turbulence-controlled, self-regulated star formation in galactic centers is expected to 
follow a scaling relation S SFR oc £ 2 . 

At very high surface densities such that t/t* ^> 1 (and assuming ^ < 1), the star 
formation rate in the radiation-pressure-dominated regime is 

^SFR,rad — 

= 710 M S kpc- yr- (g^f) " (jjj^pr) ' ( W M pc- ) ■ (W) 
which is linear rather than quadratic in the gas surface density. 

Equation ffl3|) would apply to galactic center regions experiencing all but the most 
extreme starburst activity, with a transition to the radiation-dominated regime of equation 
f lT4|) at r ~ r#, corresponding to a gas surface density 

S trans , rad = 8 x 10 3 M pc- ( - Kl f - ) " ( 1 f - e \^ d ) " . (15) 



10 cm 2 g- 1 J V3000 km s" 1 J \Q.2 x lO" 4 

The self-regulated turbulent galactic-center regime of equation ( fl3l) connects with the "mid- 
disk" regime, where internal GMC processes control star formation (see Section 1), at lower 
surface densities. Setting equation ( Fl3i) equal to S rfr = SAsf omc with tsF.GMC compa rable 
to the empirical mid-disk value tsF.moi ~ 2 Gyr (IBigiel et al.l 120081 ; iBlanc et al.l 120091 ) and 
assuming ^<1, the nominal transition to the mid-disk regime is at 

V KA A/T -2f / P*/m* \ AsF,GMc\ 1 

S trans , G MC = 54 M pc /^ 3000kmtrl j (^^j • (16) 

Physically, this transition occurs because at low surface densities, the energy input from star 
formation is insufficient to prevent "diffuse" molecular gas from becoming concentrated in 
individual gravitationally-bound GMCs. 



3. Comparison of Ssfr vs - £ to numerical simulations 

To initiate a numerical study of star formation under high surface density conditions, 
we begin with a very simple computational model, which nevertheless includes sufficient 
ingredients that we are able to investigate self-regulation of star formation in turbulence- 
dominated disks. Our simulation domain represents a local region within a starburst disk, 
which we resolve both vertically and horizontally. The domain is rotating with the disk, so 
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we include centrifugal and Coriolis forces; we assume a solid-body rotation curve, so there is 
no large-scale rotational shear within the domain. We assume that the gas cools efficiently to 
maintain an approximately constant (low) temperature, so we adopt an isothermal equation 
of state with sound speed c s = (P/p) 1//2 . For this first study, we neglect magnetic fields. 
Thus, the dynamical equations we solve are: 
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and 



V 2 $„ = AnGp 



(19) 



(e.g. iKim et al.ll2002l . taking dimensionless shear parameter q = 0). For a spherical bulge, 
the vertical gravitational field is given by g cx t = —Q 2 zz. 



We integrate the dynamical equations using a version of the Athena code (IStone et al. 



2005 



2008) , whic h employs a single-step, directionally-unsplit Godunov scheme (IGardiner fc Stone 



20081 ) to obtain conservative, second-ord er accurate solutions . The P oisson equation is 



solved using the Fourier transform method of iKoyama &: Ostrikerl (j2009al ). suitable for disk 
problems that are treated as having open boundary conditions in the vertical (z) direction 
and periodic boundary conditions in the horizontal directions. The equations above can 
be solved either for a three-dimensional dom ain, or a two-dim e nsional domain representin g 
a radial- vertical slice through the disk (e.g. IKim et al.l l2006t IKoyama fc Ostrikerl l2009al ). 
Two-dimensional radial-vertical models, which we adopt here, capture the physics of the 
competition between turbulent driving and gravitational settling in the vertical direction, 
and allow much more extensive initial exploration of parameters than is possible for three- 
dimensional models (which we intend to pursue next). 

Motivated by observations show i ng a linear relation b etween very dense gas and star 



formation (e.g. iGao fe Solomon! 120041 : IWu et al.ll2005L |2010| ). we assume that star formation 
takes place in gas above a density threshold p t hr = P n thr at a fixed efficiency per free-fall time 
£ff ( ^thr); we also assume that the total mass in stars formed per massive star is m*, and that 
each massive star injects a momentum into the ISM. Thus, when n > nthr in a given zone 
containing total mass m, the probability of a massive star formation event occurring within 
a given time step At is set to be P = AtEff(n t h I )m/[m*tff(p)}. Here, tg(p) = [3n /(32Gp)] 1 ^ 2 
is the gravitational free-fall time at density p. For every zone where a massive star formation 
event occurs, we inject a total momentum into the gas in a spherical region of radius R s h 
surrounding the event. For every feedback event associated with a massive star, the tally of 
the total mass in stars formed is increased by m*, although we do not remove any gas from 
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Fig. 2. — Sample evolution of density structure (logarithmic color scale) in the radial- vertical 
plane, for a numerical model with E = 200 M Q pc -2 . Snapshots are separated by 5 Myr. 
The domain size is L x x L z = 60 x 120 pc 2 . 

the grid (in order to maintain constant E). The computed surface star formation rate EgpR 
is the mass in stars formed per unit time within the domain, divided by the horizontal area. 

A description of the detailed numerical methods and the results of a full parameter 
survey (varying E, Q, c s , £fi(n thr ), p*/m*, R sh , n t h r , domain size L x , L z , and numerical 
resolution), will be presented separately (Shetty & Ostriker 2010, in preparation). Here, since 
we are interested in comparing to equation fflBl to test the idea of turbulent self-regulation, 
we focus on the dependence of the star formation rate E SFR on the gas surface density E. 
We fix m* = 100 M Q , p t = 3 x 10 5 M km s" 1 , i? sh = 5 pc, c s = 2 km s" 1 , and set the 
angular velocity to f2 = Myr _1 (E/10 3 M Q pc~ 2 ) (so that the Toomre parameter would be 
unity for a velocity dispersion of 7 km s -1 ). Note that the value of c s we use is larger than 
would be provided by thermal pressure at temperatures ~ 10 — 100 K; without magnetic 
fields (which are neglected in this study), shocks would become unphysically strong if c s is 
much lower. Turbulent energy is still far larger than thermal energy for the parameter range 
under consideration, so that the disk thickness is controlled by feedback-driven turbulence 
(see Section 5 and Shetty & Ostriker 2010, in preparation). The models are integrated for 
200 Myr; a quasi-steady state typically requires only a few tens of Myr to be established, 
for the parameter regime E ^ 100 M pc~ 2 under investigation. 

Figure ([2]) shows an example of typical density structure snapshots in the radial- vertical 
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plane, from a model with £ = 200 M pc 2 , taking n t hr = 5 x 10 3 cm 3 and Eg( nthr) = 0.01. 
This density threshold and efficiency factor ar e chosen based on observations of Milky 



Way and extragalactic gas and star for mation (jKrumholz fc Tan! 120071 ; lEvans et al.l 12009 



Lada et al. 


2010; 


Heiderman et al. 


2010 



ture is highly inhomogeneous due to the dominance of turbulence, which shocks and com- 
presses the gas. Effects of recent star formation events are seen as dense circular regions; 
these drive expanding shells, which are also evident. Although structure is irregular and 
highly dynamic, gas is preferentially concentrated toward the midplane due to gravity; this 
is where star formation events take place, in the densest gas. All of our models show similar 
structure and evolution, with repeated local collapse events and feedback-driven turbulent 
excitation. In some of our simulations with a sufficiently large radial domain, we find that the 
gas may collapse into a single condensation that the feedback is unable to redisperse, because 
the gravitational potential well is too deep. As we discuss in Section 5 and the Appendix, 
in real galaxies radiation pressure may play a role in dispersing strongly bound clouds and 
superclouds. Since radiation is not included in the current simulations, we consider only 
models that do not suffer this kind of global collapse. 

In spite of their highly dynamic behavior, our simulations show that a quasi-steady 
state is established in which the vertical weight of the gas is approximately balanced by 
turbulent pressure, which itself is driven by momentum injection associated with massive 
star formation. For the same model shown in Figure (j2j), Figure ([3]) shows the evolution of 
the turbulent momentum flux in the vertical direction p Q v 2 z and the characteristic vertical 
momentum injection rate to each side of the disk from star formation 0.25(p*/m*)£sFR, as 
well as the mean vertical weight of the gas disk 0.57rG£ 2 (l + x) (the value of x ~ 0-1 f° r 
this model); the respective time-averaged values of these quantities are 357, 209, and 305 
M Q km s _1 pc~ 2 Myr -1 . From equation (JTj), f p ~ 1.7, consistent with the expectation that 
the turbulent pressure is comparable to the momentum injection rate per unit area from star 
formation. Our other models show similar temporal behavior for the various momentum flux 
terms. 

Figure (jlj) shows the time-averaged values of E SFR measured from a set of simulations 
with E = 100,200,400,800 M pc -2 . For all models we set n thr = 5 x 10 3 cm~ 3 ; we show 
cases with £ff(n t hr) = 0.005, 0.01, and 0.05. The numerical results are consistent with the 
prediction given in equation (fTBl for momentum-controlled, self-regulated disk star forma- 
tion. We note that both turbulence and small-scale self-gravity are crucial in determining 
the star formation rate in these simulations. For all of the models, our adopted n t h r exceeds 
the midplane value of the density for a purely thermally supported slab (with the adopted 
value c s = 2 km s _1 ), so no star formation would take place without localized gravitational 
collapse. On the other hand, if we had not imposed a high threshold density and had simply 
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Fig. 3. — Comparison of characteristic input vertical momentum flux from star formation 
0.25(p*/m*)£sFR (squares), turbulent pressure in the vertical direction p v^ (triangles), and 
vertical weight of the gas disk 0.57rG£ 2 (l + \) (plusses), for the same model as shown in 
Fig. 02]). Each quantity is computed averaging over temporal bins of 20 Myr. 

assumed a star formation rate equal to e^T,/ tff(^midpiane) for our adopted c s and eg = 0.01 
(see Section 5), we would have obtained the same scaling with surface density Ssfr oc £ 2 as 
found with our simulations (and as predicted by equation [13]), but a much larger coefficient: 
Ssfr = £ff4GE 2 /(v^3 c s ). For the value c s = 2 km s -1 used in the present simulations, this 
would yield a factor of ~ 6 larger Esfr than equation ( 113"]) ; lower c s would yield even larger 
Ssfr- Thus, a simulation that did not include turbulent feedback but allowed gas to cool 
to low temperature (small c s ) could significantly overestimate the star formation rate. Tur- 
bulence driven by star formation keeps the volume-averaged midplane value of the density 
lower than that of a thermally-supported disk, but self-gravitating collapse of shocked, highly 
overdense local regions still occurs, leading to further (self-regulated) star formation. 
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Fig. 4. — Results of mean measured star formation rate Esfr from a set of simulations 
with a range of gas surface density S, in comparison with the prediction given in equation 
f fT3|) (solid line). For the numerical results, triangles have es(n thr ) = 0.005, squares have 
eff(^thr) = 0.01, and circles have es(n thr ) = 0.05. 



4. Comparison of Esfr vs. S to observations 

The dependence o f Ssfr on £ predicte d by the simple model of Section 2 can be com- 
pared to observations. iGenzel et al.l ( j2010af ) have recently investigated the relationship be- 
tween molecular gas and star formation in a large sample compiled from both local-Universe 
and high-redshift galaxies. The compilation includes local ULIRGs as well as normal and 
merging galaxies at z ~ 1 — 3. Gaseous surface densities in this (and other) data sets are 
based on observations of CO emission, so it is necessary to convert the integrated line emis- 
sion, I C o, to a value of E. This conversion is accomplished using a factor a = Em/Ax) 
or Xco = Nn2/Ico, where Sh2 = S/1.4 is the surface density of molecular hydrogen in 
M pc~ 2 (the factor 1.4 accounts for helium in the total surface density), Nu2 is the 
molecular hydrogen column density in cm -2 , and Ico is in units of K km s _1 so that 
-^co — 6.3 x 10 19 a in units cm _2 /( K km s _1 ) if a has units M Q pc~ 2 /( K km s _1 ). In 
Milky Way and Local Group GMCs, values of a ~ 2 — 5 (for the CO J=l-0 transition) have 
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Fig. 5. — Solution given by equation ( TIB"]) for star forma tion rate Ssfr vs. gas surface density 
X (solid line), in comparison with data compiled in iGenzel et al.l (l2010al ). In the lower 
panel, S mo i/[ M pc~ 2 ] = 1.4aJ co /[ K km s _1 pc~ 2 ] is computed assuming^ = 3.2 for 



star-for ming galaxies (stars) and a = 1 for merger systems (triangles), following IGenzel et al. 



(]2010al ). In the upper panel, a variable value of a is adopted such that a = 3.2 for £ = 
100 M pc~ 2 and a = 1 for £ = 1000 M pc -2 . Dashed lines show the least-squares fit to 
th e data, with slopes 1.4 in the lower panel and 1.9 in the upper panel. The value obtained 
by lYusef-Zadeh et all (120091 ) for the Milky Way center is also shown (squares in both panels). 



-18 - 



been inferred based on a variety of empirical techniques (e.g. IBlitz et all 120071 ). whereas 
there is evid ence that a is significantly low er in ULIRGs and other extreme star-forming 



systems (e.g. ISolomon &: Vanden Boutll2005l ). 



For surface densities £ £ 100 M© pc -2 in galaxies, CO emission can be understood as 
arising primarily from collections of gravitationally bound GMCs with limit ed variation in 
indiv i dual properties such as total column d ensity and kinetic temperature (ISolomon et al. 
19871 ; iBolatto et al.ll2008l ; iHeyer et al.ll2009l ). such that there is an approximately constant 
conversion facto r a. The standard va l ue of Xco = 2 x 10 20 based on CO emission in 
the Milky Way f lstrong fc MattoxHl996l : bame et alJl200l( l corresponds to a = 3.2. Shetty 
et al (2010) have applied radiative transfer to numerical simulations of turbulent clouds, 
confirming that mean ratio of N^/Ico is comparable to standard empirically-estimated 
values Xco = 2 x 10 20 , provided that Ay is neither very low nor very high. 

For galactic center regions with high surface densities, however, the molecular gas may 
respond to the overall gravitational potential in the disk rather than being bound in indi- 
vidual GMCs (raising the velocity dispersion), and may be systematically warmer for higher 
SgFR- These effects would incre ase the integrated C O emission J co for a given E, implying 
a smaller Xco or ol is needed ( ISolomon et al.lll997l ). Since the dependence of Esfr on S 
is nonlinear, gas surface densities may also be overestimated in some cases if the gas is not 
spatially resolved and the characteristic radius of the IR-emitting starburst is adopted as a 
surrogate (this spati al mismatch is known to be an issue for high-redshift galaxies; see e.g. 
Bothwell et al.ll2009l ). In addition, lack of spatial resolution combined with galactic inclina- 
tion can broaden observed CO line widths (due to rotation and to radial and/or azimuthal 
streaming within a beam) beyond the true turbulent velocity width in a local patch of the 
disk; without detailed modeling, it is not clear how this would affect the conversion between 
Jco an d AW For ULIRGs with typical £ ~ 10 3 M pc -2 , empirically-estimated values are 



a £ 1 



(IDownes fc Solomon! Il998l ). 



From a theoretical point of view, the value of a in a self-regulated starburst system 
should depend primarily on the disk's gas surface density S, since this sets Esfr and thus 
the gas velocity dispersion and the distributions of temperature and density through star 
formation feedback. Using numerical simulations for a range of disk properties combined with 
detailed radiative transfer calculations, it will ultimately be possible to obtain predictions for 
the dependence of a on local conditions. Empirically, a sys tematic trend for a to decrease 
with increasing S is already evident (e.g. iTacconi et al.ll2008l ). The simplest possible relation 
would be for J C o, and therefore a, to depend on £ as a power law. For example, if a is 
fixed to 1 at S = 10 3 M pc 2 from ULIRG observations, and to 3.2 at E = 10 2 M pc 2 
from the limit of "normal" GMCs, a power-law assumption yields a = (E/1000) -0 ' 52 = 
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In Figure (jSj), we compare the predicti on of equation ([TBI with the observed relation 
between E and Esfr from the compilation of iGenzel et al.l ( )2010al ). Also s hown is the surface 



densit y and star formation rate in the Milky Way Galactic center, from lYusef-Zadeh et al. 



( 120091 ). For the theoretical result, we have set x ~ ► in equation ( !T3|) from equation (JSJ), 
since C < 0.1 in the observed sample (except for two cases which have C < 0.5) if w e assum e 
the vertical velocity dispersion is ~ 10 km s _1 . In the presentation of IGenzel et al.l (j2010al ). 
a value a = 1 is adopted for merger galaxies (LIRGs and ULIRGs), while a value a = 3.2 
is adopted for other systems. In Figure (jSJ), we compare to observations using both these 
standard conversion factors (lower panel), and using the simple power-law function for a 
given above (upper panel). 

Evidently, the theoretical prediction of equation f[T3"]) follows the data over the full range 
of surface densities if a E-dependent a is used. If a = 1 is used for mergers and a = 3.2 for 
normal galaxies, equation ({TBI) agrees with the data at low E, but at high E there is an offset 
between the two clas ses o f galaxies, with only the mergers close to the result of equation (JTB"|) . 
Genzel et al.l (j2010al ) and lDaddi et al.l (120101 ) previously noted this offset, and suggested that 
a difference in dynamical timescale may be responsible for the apparent difference in 
star formation laws. With a depending on E, star formation follows a single relation for both 
mergers and normal galaxies in the E — Esfr plane, but mergers are preferentially at higher E 
- which is correlated empirically with higher Q. Regardless of the value of a, star formation 
is observed to be more efficient in merger systems, based on the ratio Li R /L co . Given the 
characteristically larger values of Jco an d S for mergers compared to normal galaxies, this is 
consistent with theoretical expectations that Esfr depends nonlinearly on E - in particular, 
as E SF R oc E 2 based on equation (jTBl . 



5. Disk properties in the turbulence-dominated regime 

In Section 2, we showed that if turbulence driven by star formation feedback dominates 
other forms of pressure, the equilibrium star formation rate is expected to be proportional 
to the weie ht of the ISMH The force balance (or momentum conservation) argument leads 



3 The star formation rate is also expected to be proportional to the weight of the ISM if thermal pressure 
domi nates other forms of pressure, in a medium dominated by diffuse atomic gas and heated by young 
stars I Ostriker et al. 201C)I ). If, however, radiation pressure dominates both turbulent and thermal pressure, 
the star formation rate would be proportional to the vertical gravitational field rather than the weight (i.e. 
reduced by a factor cx 1/r oc 1/E - cf. equations [13] and [14]). 
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to a relation Esfr oc S 2 for gas-dominated regions in a self-regulated steady state, but it 
says nothing about the internal processes or character of the disk. In order to maintain 
equilibrium, the disk must internally adjust itself so that stars are produced at the required 
rate, and this internal arrangement must also be consistent with the detailed dynamics of a 
self-gravitating, turbulent hydrodynamic (or hydromagnetic) flow. 

In a region of mean gas density p (and where p* <C p), the gravitational free-fall 
time at the mean density, t s (p) = [37r/(32Gp)] 1 / 2 , represents a characteristic dynamical 
time. The star formation rate can then be expressed as M* = £ff(p)M gas / iff(p), where 
£ff(p) is an efficiency factor that depends on the detailed properties of the gas. Although 
it is uncertain exactly what controls the rate of localized gas collapse in highly-tu r bulen t 



systems, evidence from current observations (Krumho 



theory (IKrumholz &: 



( jPadoan fc Nordlund 



2009; 



VlcKe jlgooi IPadoan fe Norxllund 



z fc Tad 120071 lEvans et all 120091 ). 
2009j ). and numerical simulations 



iRosas- Guevara et al.l I2OO9I ; Ivazquez-Semadeni et al.l I2OO9I ) sug 



gests that the star formation efficiency facto r is low, e$[p) ~ - 01> fo r molecular gas in 
observed star-forming regions at a range of p. IKrumholz fc McKed ( 120051 ) proposed that for 
turbulence-dominated systems, £g(p) is given by the fraction of mass compressed to densities 
higher than the mean by a factor ~ M 2 for M. = tWbA'th the Mach number. Because th e 
density PDF in turbulent regions is a lognormal function of Ai (IMcKee & Ostrikerl 120071 ) . 
Krumholz fc McKed (120051 ) show that the resulting dependence of the star-forming efficiency 
on the Mach number is weak, eg(p) oc A4~ ' 3 . Thus, if localized collapse to make stars is 
controlled primarily by turbulence, the overall rate is expected to be governed primarily by 
the free-fall time on large scales. 

In (self-gravitating) disks, it is convenient to use the mean midplane density po = pno 
as a reference value, so that the star formation rate per unit area can be written as 

E 



J SFR 



£ff(Po) 



ts{po 



(20) 



For a Gaussian vertical density distribution, p = po/v^ so that eg(p ) = 0.8eg(p). As in 
Section 2, in adopting the form ( 1201) . we have implicitly assu med that the gas collapses only 
locally -on scales small compared to the Toomre wavelength ( Toomrelll964 ) At = GE(27r / k) 2 
- so that angular momentum and shear do not dominate the dynamics of star formation. In 
addition, equation ( 120|) implicitly assumes that the free-fall time at the mean midplane gas 
density is shorter than the vertical oscillation time in the gravitational field of the stellar 
bulge; this will be true if gas self-gravity dominates vertical potential gradients. 

We note that equation ( 120]) is not expected to apply in regions of a galactic disk that 
have a fundamentally multiphase character, in terms of the mass- and volume- fractions of 
gas at different temperatures. In particular, if most of the ISM's mass is in clouds that 
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are self-gravitating entities isolated from their surroundings, the local density within these 
self-gravitating, star-forming clouds need not be proportional to the (volume-weighted) mean 
midplane density of the disk. Using multiphase disk simulations, iKoyama &: Ostrikerl (j2009af ) 
indeed found that equation f[2"01 would yield a different (steeper) star formation law than is 
obtained by adopting a constant star formation timescale only in very dense (gravitationally- 
bound) gas. 



For outer-disk regions dominated in mass by diffuse atomic gas, lOstriker et all (120101 ) 
argued that a key factor regulating star formation is the heating and cooling of the warm, 
volume-filling medium. Gas that cannot be maintained in thermal equilibrium in the diffuse 
phase cools and drops out to make gravitationally-bound, star-forming clouds. This increases 
the heating rate (from stellar UV), and decreases the cooling rate (since the density is 
lower), for the remaining diffuse gas. Mass dropout continues until heating balances cooling. 
In a state of simultaneous thermal and dynamical equilibrium, the star formation rate for 
an outer-disk region will still vary inversely with the dynamical time of the dominant gas 
component, analogous to equation (1201) . but the relevant dynamical time is not ts(po). In 
outer-disk regions, the ISM scale height is generally set by the gravity of stars plus dark 
matter, rather than the self-gravity of the gas. As a consequence, the gas disk's internal 
dynamical time tdyn ~ H/a z (or the timescale for vertical oscillations ~ 2itH/<j z ) varies as 
^dyn oc (Gp*) _1//2 rather than tdyn oc ts(po) oc (G^p) " 1 / 2 , yielding an outer -disk star formation 
rate S SFR oc Hy/p^ (see equations [22] and [A16] in lOstriker et al.l ( |2010| )). 



For turbulent, molecular-gas-dominated starburst regions where equation f[2"U|) is hy- 
pothesized to apply, the value of es(po) is uncertain, but analysis of observations including 
nearby star-forming clouds, and both moderate-density a nd high-density extragalactic molec 
ular gas, is consistent with a range e«{po) ~ 0.003—0.06 (IKrumholz fc Tanll2007l ; iBigiel et al. 



2008 



Evans et al. 2009; Blanc et al. 



20091 ) . Simulations of turbulent, star-forming disks can 



also be used to evaluate £e(po)- The mean midplane density p can be measured and used 
to obtain a mean value of tff(po); if the value of Esfr measured in the simulations is mul- 
tiplied by the measured tg(po) and divided by S, the result is a numerical evaluation of 
£ff(p ). For the numerical simulations described in Section 3, the models yield a range 
£ff(Po) ~ 0.004 — 0.01 (for further discussion, see Shetty & Ostriker 2010, in preparation). 
As a fiducial value, we shall adopt eg(p ) = 0.005, but it is important to keep in mind that 
further observations and simulations are needed to determine this value (and any dependence 
on parameters) more conclusively. 

A disk that is in vertical equilibrium will obey equation ([3]); assuming gas thermal+turbulent 
pressure exceeds other forms of support (7Z <C 1), po = 7rGS 2 (l + x)/(2a 2 ). Substituting 
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this expression for po in equation f l20|) . the star formation rate can be written as: 
45 ff (p )(l + x) 1/2 G£ 2 



J SFR 



^3 Vz 
0.051 M kpc~ 2 yr _1 (l + x) 1 ' 2 



gff(po 

0.005 7 VlOkms- 1 ^ V 100 M P c_2 
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Thus, disks dominated by self- rather than external-gravity < 1) will have £sfr oc £ 2 
if £fr(po) ~ constant and the velocity dispersion varies only weakly with £ - as we verify 
below. 

We now focus on the case in which turbulence driven by star formation dominates the 
vertical support and control of star formation. In this limit, equating (TT3"]) with (I2ip for 
o"z — > v z yields: 

2/pgff(po)p* 1 
Vz ~ VSnm* (l + x) 1/2 

-i f P ( £s(po)\ ( P*/m* 



= 5 - 5kms (22) 

The only dependence of the velocity dispersion on the mean gaseous surface density £ is 
through x'i if the disk is marginally unstable (Q £ 2), however, x ~ 0-5 (see Section 2). 
In the supernova-driven, turbulence-dominated limit, the velocity dispersion is therefore 
expected to vary only weakly with £ for disks that are marginally gravitationally unstable. 
The numerical simulations described in Section 3 indeed show a very weak dependence of 
v z on £, with turbulent velocity dispersions varying by a factor two (~ 4 — 9 km s _1 ) in 
models with gas surface density between £ = 100 M pc -2 and £ = 800 M pc -2 , which 
have £sfr varying by more than a factor 100. The main parameter controlling the velocity 
dispersion, in disks that cool strongly and therefore have highly-supersonic turbulence, is the 
specific momentum input rate associated with star formation, p*/m*. In equation (122]) . the 
fiducial value of p*/m* chosen is that associated with cooled supernova shells. In addition 
to the direct momentum input from radiative blast waves, cosmic rays that are accelerated 
by blast waves may drive additional turbulence in escaping from the disk (e.g. via Parker 
instabilities), which would raise the value of and increase v z . 

By combining equations ( 122"]) and (JS]) (for TZ <C 1), we can solve for x and C = x(l + x) 
in terms of the star formation parameters and the gas disk and bulge properties. Defining 

2Q ( 4/p£ ff (p )p, y p b 

7r y 37rm, J G£ 2 ^ 
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= 011 ( E « M Y( p ' /m ' Y( " Y( ? V 2 (24) 

V 0.005 J VSOOOkms- 1 ; \0.1 Myr" 1 / \100 M pr^ v ; 

we obtain the implicit relation x(l + x) 2 — <A. An approximate solution good to within 15% 
over all A is 

X - T7~p73- (25) 

For .4. <C 1, x ~ C ~ «4j an d t ne Toomre parameter is Q = 2.5(aji/a z )A 1 ^ 2 . 

In the turbulence-dominated limit o z rs v z (and for 72. <C 1), the mean midplane density 
p = no(1.4m#) (for m# the hydrogen mass) and half-thickness 77 = X/(2p ) are given by: 

(1 + x)7r GS 2 3tt 2 / im» \ 2 2 
Po = n — = ^r( 1 + X) 757 GS 26 



so that 



-si 1 + x)V £ ff(^o)\ 2 ( P*/m* x 



ri0 = 64cm_ — ^IcWj U000 kms-i j (, 100 M pc -' ' ' (27) 



and 

77 



2 1 1 /or c \ \ 2 



1 + x) ttCS 2 3vr (1 + x) 2 V 7rm» / GE 

23 pc fp f ^(po) \ 2 f P*/m* V f ^ | K 

1 ;i + x ) 2 V 0.005 7 VSOOOkms-V V!00 M pc- 2 ' 

We have written equations (TI3]) . ( )22l) . ( J26l) . and (J28l) in such a way as to highlight the de- 
pendence on the gas surface density. As £ increases, and taking A, x "C 1, the surface density 
of star formation increases as S SFR oc E 2 m*/p^,, the vertical velocity dispersion is constant 
t>2 oc £fr(po)p*/ m *, the mean midplane gas density increases as po oc E 2 ( £ff(po)p*/ m *) -2 and 
the scale height decreases as 77 oc S _1 ( eg (po)p*/m*) 2 . Increasing the specific momentum 
input Px/m* from star formation decreases the star formation rate and increases the veloc- 
ity dispersion and disk thickness. The absence of any explicit timescale in equation (fT3"l) 
is a signature of self-regulation: the star formation rate must adjust until the turbulence 
driven by feedback from young stars provides a pressure that matches the vertical weight 
of the gas. If the turbulence level were lower than the equilibrium value, the mean density 
would be higher than equilibrium (p oc v~ 2 ), which would then lead to a shorter dynamical 
time tjf oc p -1 / 2 and consequently (from equations [20] and [2T] ) a higher star formation 
rate, driving the velocity dispersion upwards. Similarly, too high a turbulence level would 
lead to lower-than-equilibrium p and Xsfr, which would reduce turbulent driving from star 
formation and eventually lead to a lower value of v z . 



-24- 



Although a balance between turbulent driving and dissipation has not been explicitly 
used in order to derive the above results, it is straightforward to see that these consider- 
ations yield equivalent results. The turbulent vertical momentum per unit time per unit 
gas mass that is driven by stellar inputs is ~ p t E SFR /(Sm l! ). Since the d issipation time for 



turbulence is compar able to the flow crossing time over the largest scale (IStone et al.l Il998 



Mac Low et al.lll998l ). which in this case is the disk thickness H ~ v 2 /(GH), the dissipation 
rate of momentum per unit gas mass is ~ v 2 z /H ~ GS. Equating driving with dissipation 
yields S SFR ~ GT^m^/p^, which is the same as equation (|T5|) up to order-unity dimensionless 
constants and the factor f p that has been incorporated to parameterize the details of the 
momentum injection and mixing. 

In terms of the feedback and disk parameters, the Toomre parameter for turbulence- 
dominated disks is 

V ' (l + X y /2 v z V 0-005 J \3000kms-V \0.1 Myr"V K 1 ^ M © P^ 2 / ' 

Thus, for £ and Q in the range observed for galactic center starbursts, the disk will be near 
the margin of instability. Even if Q ^ 2, the intermittency of turbulence implies that some 
regions may temporarily become unstable. If a large region of the disk becomes unstable 
and begins to contract, a gravitationally bound supercloud may form. 

In a gravitationally bound supercloud, the outward force F^ n from direct kinetic feed- 
back from star formation varies as ~ M*p*/m*, whereas the outward force F rac j from repro- 
cessed radiation varies as ~ L*r/c oc M*M/r 2 , where M#, L*, M, and r are the supercloud's 
star formation rate, luminosity, mass, and radius. The inward gravitational force -Fg rav varies 
as ~ GM 2 /r 2 . If M* oc Mj tg(p) oc r~ 3 / 2 , then Fki n /F grav oc r 1//2 (decreasing as a supercloud 
shrinks), whereas -F ra d/-F grav oc r -3 / 2 (increasing as a supercloud shrinks). This suggests that 
superclouds, if they become bound, could be destroyed by radiation pressure feedback but 
not by kinetic feedback. In sufficiently optically thick disks, radiation pressure might prevent 
supercloud formation; a simple estimate in Appendix B suggests that disks with gas surface 
densities up to S ~ (3000/kir) M pc -2 may be susceptible to supercloud formation. 



6. Summary and discussion 

Increasingly high-resolution, multiwavelength observations of galaxies have refined em- 
pirical relations between the gas content and mean star formation rate in galaxies. The 
picture now emerging appears to include at least three regimes: outer-disks, mid-disks, and 
galactic centers. The first and last regimes evidence superlinear mean relations between the 



-25 - 



surface density of star formation Esfr and the gas surface density E, while the relation is 
close to linear for mid-disks. 

From a theoretical point of view, star for mation in all regimes is likely to be self- 



regulated in some way. For the outer-disk regime, lOstriker et al.l ( 120101 ) have proposed that 
UV feedback from massive young stars is particularly important to self-regulation. In this 
model, the star formation rate in outer disks must be such that the heating of diffuse atomic 
gas by stellar UV radiation balances cooling, with the cooling rate set by the gas pressure (and 
hence responsive to the vertical gravity of the disk). For the mid-disk regime, gas resides 
predominantly in gravitationally-bound (but transient) GMCs with mean densities (and 
hence star formation rates) that depend on the amplitude of supersonic turbulence driven 
by star formation feedback (possibly by expanding HII regions, although this is not well 
understood). For galactic center regions, we argue in this paper that the momentum injected 
to the disk by massive star formation is crucial to establishing a self-regulated state. Vertical 
equilibrium of the disk requires gravity to be balanced by an outward momentum flux, 
consisting primarily of turbulent pressure (largely driven by supernovae) for most molecule- 
dominated central regions. Our analysis suggests that the galactic-center star formation 
regime may further subdivide, with a transition to radiation-dominated momentum flux (or 
pressure) only for extremely high E, highly-opaque disks. 

The analysis of Section 2 provides a prediction (see equation [10]) for the dependence 
of EgpR on E and the density p& of the stellar bulge, parameterized by the mean specific 
momentum p^/m* injected by star formation to the ISM (primarily from supernovae), by 
the mean IR opacity Kjr (which determines how much reprocessed starlight is trapped), 
and by the mass-to-radiation energy efficiency of star formation e* (which depends on the 
IMF). Assuming that gas gravity dominates the vertical potential, self-regulated, turbulence- 
dominated disks are expected to have E SFR ps 27rm* GE 2 /^ (see equation [TS]), whereas 
radiation-dominated disks are expected to have S SFR rs 27rGE/(e*CKrR) (see equation |T3|) . 
The turbulence-dominated regime is expected to apply for E ~ 100 — 10 4 M Q pc -2 , covering 
most starbursts (although not the innermost disk regions that merge into AGN disks). For 
the radiation-dominated limit, an equivalent expression to our result was previously obtained 



(under somewhat different assumptions) by iThompson et al.l ( 120051 ) (see their equation 28 



and they also noted that optically-thin disks are expected to have a steeper dependence of 
Esfr on E than optically-thick disks. 

In a self-regulated disk, the feedback parameters determine the star formation rate, 
independent of the internal structure of the disk. Thus, the vertically-integrated star forma- 
tion rates of equations (|T3|) and (JHJ) depend only on the vertically-integrated disk prop- 
erties. However, it is natural that the star formation rate should also connect to the 
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timescale for gas to become concentrated into high-density clouds, as discussed in Section 
5. For a system in equilibrium where the gas is confined primarily by its own gravita- 
tional potential, this timescale is the gravitational free-fall time at the mean density; for 
a disk, 1/ tff(p) oc (GE/if) 1 / 2 . The mean density (or disk thickness), in turn, depends on 
what is balancing gravity - thermal pressure, turbulent pressure, or radiation pressure. In 
the turbulence-dominated, gas-gravity dominated case, H = t> 2 /(7rGE) for v z the vertical 
velocity dispersion, which yields l/tg(p) = 4GE/(v / 3f 2 )- If £s(Po) is the collapse effi- 
ciency per free-fall time at the mean midplane density, the star formation rate is given by 
^sfr = 2.3 es(po)GT l 2 /v z (see equation l2Tj) . 



Krumholz fc McKed (120051 ) proposed that the log-normal density distribution in highly 
turbulent systems will lead to a star-forming efficiency declining as the —1/3 power of the 
turbulent Mach number. Including this scaling would lead to Esf r oc E 2 tC 1,3 . Under the 
assumption that in molecule-dominated gas Esfr = £ff(p)S/ t$(p), IKrumholz et al.l (120091 ) 
proposed a power-law relation between Esfr anc ^ ^ ^ n galactic center regions. In obtaining 
this relation, they assumed that the Toomre parameter Q f=a 1, and also that Q oc E 0,5 ; 
together these would imply that velocity dispersion varies oc E a5 . Although we in fact 
argue that the vertical velocity dispersion v z is rather insensitive to E, inser ting v z oc E 0,5 ^ 
Erfr oc E 2 u~ 1,3 would yield Esfr oc E 13 , the result given in equation (10) of IKrumholz et al. 

(boool). 



In order for the disk to form stars at a rate controlled by the free-fall time (Esfr ~ 
2.3 £ff(po)GE 2 /t> 2 ), and also for the momentum feedback from star formation to control 
vertical equilibrium of the disk (Esfr ~ 27r m *GE 2 / p*) , the vertical velocity dispersion must 
adjust to v z ~ OAEff(po)p*/m* (see equation [22]). Interestingly, this is independent of the disk 
surface density and the star formation rate; our numerical simulations in fact show that v z is 
insensitive to the star formation rate. Other recent numerical simulations have also indicated 



(Dib et al. 


2006; 


Shettv & Ostriker 


2008; 


Aeertz et al. 


2009; 


Joune et al. 


2009) 



v z ~ ^ff(po)p*/^* can also be understood as simply a balance between driving turbulent 
velocities in the ISM at a rate ~ £fr(po)/[ m * ^ff(Po)]; an d dissipating turbulent velocities 
at a rate ~ v z v z /H ~ v z j t s (p ). For es(p ) ~ 0.01 and p*/m* ~ 3000 km s _1 , the velocity 
dispersion driven by star formation feedback is v z ~ 10 km s _1 . 

More generally, if the vertical momentum injection rate per unit ISM gas mass from 
star formation is ~ (1 — 2)p*/(4m !| ,£sF,gas), and the corresponding turbulent dissipation rate 
is ~ v z /[(l — 2)t dyn ] for vertical dynamical time td yn = H/a z , the expected vertical velocity 
dispersion is v z ~ (0.3 — l)(p*/m*)td y n/^SF,gas- This is within a factor of a few of 10 km s -1 if 
p*/m* ~ 3000 km s^ 1 from supernovae and tsF, g as ~ lOOtdyn (which is true in outer disks as 
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well as starbursts; see below), regardless of whether the vertical dynamical time is set by gas 
self-gravity or by gravity of an external potential (stellar bulge, stellar disk, or dark matter 
halo). Velocity dis persions reported for starburst regions are of ten much larger than the 
values we find (e.g. iDownes Solomonll 19981 ; iGenzel et al.ll2010bl ). but it is difficult to elim- 



inate sub-beam sheared rotation, radial/azimuthal streaming, and other non-turbulent (and 
non- vertical) motions on the relevant scales (~ H <C 100 pc). Observations of molecular 
yelocity dispersions within t he central kpc of a few nearby, face-on galaxies are ~ 10 km s _1 
( Combes fc Becauaert 1997T); tu rbulent molecular velocity dispersions from larger scales are 



also similar ( jWilson et al.l 120101 ) . 



The momentum-regulated, turbulence-dominated prediction of equation ffTB"]) for Esfr 
can be compared to the results of numerical simulations and to observations. We do this in 
Sections 3 and 4, respectively. The simulations indeed yield a relation consistent with the 
prediction, £sfr ~ 0.1 M kpc~ 2 yr _1 (S/100 M pc~ 2 ) 2 , adopting a value of the momen- 
tum input per unit stellar mass formed of p*/m* = 3000 km s _1 . Velocity dispersions from 
the simulations are also comparable to prediction. Further comparisons of the star formation 
rates and disk properties from the model with the results of numerical simulations will be 
presented in Shetty & Ostriker (2010, in preparation). 

A difficulty in comparing to observations is that the molecular mass must be obtained 
indirectly through observations of CO transitions. Although the conversion factor Aco from 
integrated CO intensity to hydrogen column density has been well calibrated empirically 
using several different methods for individual GMCs and main-disk regions which integrate 
over collections of GMCs, the value of Aco is much less certain in galactic center regions. 
There is evidence that Aco decreases in high-S starburst regions; this may owe to a combi- 
nation of factors, including higher CO excitation in the warmer, denser gas, a larger fraction 
of the gas that is sufficiently shielded to create CO, and larger velocity gradients that allow 
ra diation to escape mo re easily. We compare to the large sample of observations compiled 



in 



Genzel et al.l ( I2010af ). If the conversion factor is assumed to decrease with increasing in- 
tegrated CO intensity (as Aco oc /co' 3 ), ^ e observed and predicted relations for Ssfr vs. S 
are in excellent agreement, for the range E ~ 10 2 — 10 4 M pc -2 . If two different constant 
conversion factors are adopted for normal and merger systems, the former agree with the 
predicted Ssfr at low E, and the latter agree at high S. To make progress in relating theory 
to observations, a top priority is to determine how Xco varies under the conditions that 
prevail in galactic center regions. This will require both empirical calibrations, and radia- 
tive transfer models applied to numerical simulations that resolve the structure of turbulent 
disks. 



Finally, we note that the star formation timescale tsF, g as = S/Ssfr ~ p*/(27rGm*£) in 
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self-regulated, turbulence-dominated, molecular galactic center disks can als o be expressed 
as tsF. oas ~ ^dyn/ £ b1po) ~ lOOtdyn, for vertical dynamical time t dyn = H/a z . In lOstriker et al. 



( 12010| ). we showeccl that in outer disks, self- regulation via UV heating of atomic gas leads to 
a quantitatively similar result, tsF.gas ~ 200H/a z . The constant factor ~ 200 in the outer- 
disk formula for tsF, g as is set, however, by the ratio of heating and cooling rate coefficients 
for atomic gas and by the UV production efficiency of the young-star population, rather 
than by the turbulent processes that set the constant factor ~ 1/ £ff(po) i n the galactic- 
center formula. Thus, although very different physical processes control star formation in 
galactic centers and outer disks, the dependence of Esfr on 1/tdyn (and on Q ~ Q/tdyn) is 
similar, both in terms of the scaling and (coincidentally) the coefficient. We believe that this 
dependence on 1/tdyn, together with the limited range of Q that disk evolution yields, is the 
reason that both whole galaxies and cent ral starbursts, at a range of redshifts, t end to follow 
the em pirical relation Esfr ~ 0.02£f2 ( IKennicuttl Il998bt iDaddi et al.l |2010| ; iGenzel et al. 



2010af ). 
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A. Appendix: momentum injection to the ISM from radiation-driven shells 

around star clusters 

Consider a GMC of mass Mqmc which forms stars with a net efficiency over its lifetime 
of £gmc = M*/M GMC , ejecting the remaining gas mass (1 - £gmc)^gmc- Treating the 
stellar component as collected in a single central cluster of (fixed) radius r c i and luminosity 
L* = Sgmc^A^gmcj and the ejected mass as an expanding spherical shell of (variable) 
radius r surrounding it, the total radiation force applied to the interior of the shell by 
reprocessed (diffuse) radiation is F rad = 47rr 2 (a ra d/3)T i 4 lt , where T^ t = (3/4)rT e 4 ff for shell 
optical depth r ^> 1 and internal and surface (effective) temperatures T int and T eff . Equating 
L* with 4vrr 2 (a rad c/4)r e 4 ff and using r = Ki R (l-£:GMc)^GMc/(4vrr 2 ), F rad = L^r/c = k ir (1- 
£gmc)£gmc^Mq MC / (An cr 2 ). The gravitational force on the shell, including that of the 



See equation (A16) in lOstriker et al. ( 2010 ): note that i C on there is the same as tsF,gas here. 
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central cluster and its own self-gravity, is — GM GMC (1 — ^gmc)/^ 7 " 2 ) 
of the shell is therefore 



The net acceleration 



GM, 



GMC 



^IR^GMC (1 + £GMc) 



AtiGc 



(4>£gmc - 1) GM GMC 



(Al) 



for if) = ^Kir/ {2ncG) — 1. As pointed out by lMurray et al.l (120101 ). this implies that there is 
a minimum value of £gmc f° r the pressure associated with reprocessed radiation to disrupt 
the cloud; the shell can only become unbound if £gmc > = £mm- 



Equation (I A 1 [) can be integrated to obtain the asymptotic speed Vf = [(V^gmc — 
l)GMQMc/r c i\ 1 / 2 of the shell, assuming its initial radius is r c \. The resulting ratio of the 
momentum in the shell divided by the total mass in stars formed, representing p*/m* for 
a radiation-driven shell, is vj(l — £gmc)/ £ gmc- F° r anv given value of ip, it is straight- 
forward to show that the maximum momentum-to-mass rat io is obtained whe n £gmc — 
Ae min [l + (l+%e min ) l l 2 }- 1 = e max . Usin g ^ = 20 erg s' 1 g" 1 JDopita et al. l l2006h for young, 
luminous clu sters that fully sample thelKroupal (20011) IMF and Kir = 20 cm 2 g _1 for warm, 



dusty shells (IMurray et al.l |201Q| ; ISemenov et al. 



20031 ). this yields e r 



0.5, 



0.6, 



and [p*/m,] rad = QA{GM GMC /r d ) l l 2 = 30 km s^W/lO 6 M & ) 1/2 {r cl / pc)" 1 / 2 . This is 
comparable to the velocity dispersion of the cluster itself, which even for t he most massive 



clust e rs in starbursts does not exceed 100 km s 1 and is typically much lower (IMcCrady et al. 
200,4 iMaraston etaDbooi bverzier et aDboOol ). 



In addition to diffuse (reprocessed) radiation, streaming radiation (dominated by UV) 
from the star cluster also imparts momentum to the expanding gas cloud (or shell) where 
it is first absorbed. The maximum force from this streaming radiation is L*/c, such that 
the maximum contribution to Vf{l — £gmc)/^gmc — V*/ m * is / ^fdt/c = e*c ~ 200 km s _1 . 
Thus, allowing for both reprocessed and direct stellar radiation, the value of p*/m* from 
radiation-driven expanding shells is smaller than the momentum-per-stellar mass injected to 
the disk by radiative supernovae shells, ~ 3000 km s _1 . 



B. Appendix: supercloud evolution 

Potentially, large-scale instabilities in the disk can lead to gaseous supercloud formation. 
How would a supercloud evolve if it is able to form? For a spherical, uniform- density cloud of 
mass M and radius r that forms stars at a rate M* = es(p)M/ tg(p), producing a luminosity 
= e*c 2 M* and a supernova rate M*/m*, the ratio of the outward force (due to the 
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combined momentum input from supernovae and reprocessed radiation pressur^f]) to the 
inward force (due to gravity) is 



?>tiG 1 I 2 M 1 I 2 



f p * r.i/2 i 3e«c/ti R M 3/2 



47T 



(Bl) 



The contribution to F + /F_ from supernovae declines oc r 1//2 as r decreases, whereas the 
contribution to F + /F_ from radiation increases oc r~ 3 / 2 . Thus, radiation forces become 
increasingly important as the cloud shrinks; eventually, if F + /F_ becomes > 1, inward con- 
traction will be halted. If the star formation efficiency is sufficient, the cloud will ultimately 
be destroyed, with a substantial fraction of its mass once again becom ing unbound (cf. 
Krumholz fc Matznerl l2009t Murray et al.l I201CH : iKrumholz &: Dekell l2010l and the discussion 
in Appendix A). 

The masses of clouds tha t form via gravitational instability in disks w ith Q < Q cr it are 



typically 
(7r#) 2 £ = 



Kim et al.ll2002l . 120031 : iKim fc Ostrikerl 120071 ) . where M Ji2 d 



^ 1 - 10M Ji2D (e.g. 

cr 4 /(G 2 E) for a the gas velocity dispersion. Letting M = mMj j2 D, r = %H, and 



assuming a = v z , we can use equations 

F+ l.| 



m 



1/2 



x l ' 2 + 



, ((IIP, and (d2D to obtain 



3ir t m 



4 n x 3 / 2 



(B2) 



where r is evaluated using the unperturbed disk surface density E, and we have assumed 

For a given m, the force ratio has a minimum at x m ; n = 2.7(mr/r t ) 1 / 2 , where its value 
is (F + /F_) min = 4.0m" 1 / 4 (r/r t ) 1 / 4 . Taken at face value, this suggests that only disks with a 
sufficiently low surface density £ will form contracting superclouds, since only if r/r* is suf- 
ficiently small can F + /F„ be < 1. With m = 10, F + /F_ < 1 for t/t* < 0.04, corresponding 
to S < (3000/k ir ) M pc~ 2 . Note also that x min < 1 only for £ < (1000/k ir ) M pc" 2 . Of 
course, given the highly idealized assumptions adopted, the particular value obtained from 
this simple estimate should not be taken too seriously. Physically, however, it is reason- 
able to expect higher-E disks that are increasingly radiation-pressure supported to be less 
subject to gravitational instability, because of their stiff equation of state: if P ra( j oc L^p/r 
and oc Mp 1//2 , the internal radiation pressure in a supercloud would vary as P oc p 11 / 6 . 
Potentially, more detailed understanding of large-scale gra yitational in s tabili ties in starburst 
disks could be obtained via an analysis similar to that of iThompsonl (120081 ). but including 
sources of radiation and turbulence. 



5 Streaming radiation adds a term < e* c ~ 200 km s 1 to p^/m* that is small compared to the contribution 
Po/to* ~ 3000 km s _1 from supernovae, and thus may be neglected. 
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